Effects of Fermi energy, dot size and leads width on weak localization 

in chaotic quantum dots 
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Magnetotransport in chaotic quantum dots at low magnetic fields is investigated by means of a tight 
binding Hamiltonian on L x L clusters of the square lattice. Chaoticity is induced by introducing L 
bulk vacancies. The dependence of weak localization on the Fermi energy, dot size and leads width 
is investigated in detail and the results compared with those of previous analyses, in particular with 
random matrix theory predictions. Our results indicate that the dependence of the critical flux $ c 
on the square root of the number of open modes, as predicted by random matrix theory, is obscured 
by the strong energy dependence of the proportionality constant. Instead, the size dep endence of 
the critical flux predicted by Efetov and random matrix theory, namely, <!> c cx y^l/L, is clearly 
illustrated by the present results. Our numerical results do also show that the weak localization 
term significantly decreases as the leads width W approaches L. However, calculations for W — L 
indicate that the weak localization effect does not disappear as L increases. 
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I. INTRODUCTION 

Experimental studies of magnetoconductance in quan- 
tum dots show that, at low magnetic fields (typically be- 
low one flux quantum), the conductance increases with 
the field EHl . The effect has been investigated theoreti- 
cally d||l||| an d related to a similar behavior observed 
in disordered metallic conductors in the diffusive regime, 
that is referred to as Weak Localization (WL) P,[lO||. 

There is also fairly conclusive experimental evidence 
which indicates that the average magnetoconductance 
< G(B) > behaves in a qualitatively different way in reg- 
ular and chaotic cavities, namely, whereas in the former 
it increases linearly with B, in chaotic cavities the WL 
peak has a Lorentzian shape Semiclassical analyses 
ascribe this difference to the distributions of the areas A 
enclosed by the trajectories of the carriers §. While m 
regular systems the probability distribution of enclosed 
areas larger than A is cx I /A [|TI|, in fully chaotic sys- 
tems it is exponential As a consequence, in chaotic 
cavities the increment in the magnetoconductance as a 
function of magnetic flux $ is given by: 



SG = G($) - G(0) 



a$ 2 



I + 6$ 2 



(1) 



where the conductance and the magnetic flux are given 



in units of their respective quanta, Go 



- /h and 



$o = h/e. The constant b gives the critical flux at which 
the time-reversal symmetry is effectively destroyed, $ c = 
1/Vb, whereas the ratio a/b gives the weak localization 
term, i.e., G(oo) — G(0) = a/b. The supersymmetric 
cr-model predicts that a and b should be inversely pro- 



portional to the number of channels N c h that contribute 
to the current [flOl, 



vD L 



(2) 



where c is a constant (for the present geometry c — 2iv/3), 
v the density of states, Do the diffusion coefficient and L 
the linear size of the cavity. The size dependence arises 
from the standard expression for the diffusion coefficient 
D Q = vpl/2, where vp is the Fermi velocity, and I the 
clastic mean free path, and from the fact that in a two- 
dimensional ballistic system, I cx L p3| . The qualitative 
behavior of Eq. (||) is similar to the Random Matrix The- 
ory (RMT) result for the critical flux at which the time- 
reversal symmetry is broken (GOE-GUE transition) re- 
ported in Q . The two constants a and b are proportional 
to each other. In particular RMT gives 



(3) 



b 4N ch + 2 ' 
where N c h is related to the zero field conductance 



through, 



G 



RMT 



(0) 



Aeh 

2 



N. 



ch 



(4) 

On the other hand, a fitting of the numerical results ob- 
tained from a random matrix model Hamiltonian gave 



b = 2k- 



2N r „ 



1 



(5) 



k being a constant which, as in Eq. (Q), depends on the 
Fermi energy. Although Eq. (§) gives the same depen- 
dence on the number of channels than Eq. (j^) in the large 
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N c h limit, it does not explicitly reproduce neither its size 
nor its energy dependence. Moreover, as remarked in 
Q |; Eq. (||) is only valid for few channel ballistic cavities. 
It should also be mentioned that the size dependence of 
the constant b has also been obtained within RMT (see 

EMU). 

At present there is no published numerical study of the 
effects of the size of the cavity, the leads width, and the 
Fermi energy on weak localization in reasonably realistic 
models of quantum chaotic cavities. The purpose of this 
work is to discuss the results of such an investigation. 
Quantum dots are described by means of a tight-binding 
Hamiltonian on L x L clusters of the square lattice. Non- 
regular (chaotic) behavior is induced by introducing a 
number of bulk vacancies proportional to the linear size 
of the system fl7]| . This model has been shown to be- 
have similarly to dots in which chaoticity is induced by 
introducing disorder at the surface [ ^8|Jl9[| . The effects 
of leads width W, system size, and number of channels 
that contribute to the current are discussed in detail. 
Our results show that the the critical flux is not simply 
proportional to the square root of the number of open 
channels as concluded in juj; it turns out that this rela- 
tionship is obscured by the strong energy dependence of 
the proportionality constant already implicit in Eq. (|^). 
Significant deviations from RMT are observed for large 
leads width (W of the order of the system size L). In 
particular the weak localization term decreases as W ap- 
proaches L. However, our numerical data for W = L 
indicate that this term does not vanish as L increases. 

The paper is organized as follows. Section II includes 
a description of our model of chaotic quantum dot and 
of the method we used to compute the current. The re- 
sults are discussed in Section III. We first briefly consider 
the case of zero field, comparing our results with those 
derived from RMT. The results concerning the effects 
of Fermi energy, leads width and dot size are presented 
and discussed thereafter. Again, comparison with RMT 
is highlighted. Section IV is devoted to summarize the 
conclusions of our work. 



II. MODEL AND PROCEDURES 

A. Model of quantum chaotic dot 

Our model of a quantum chaotic dot is described 
by means of a tight-binding Hamiltonian with a single 
atomic orbital per lattice site, 



H 



j m,n:m' ,n 



|m, n >< m', n |, 



(6) 



where \m, n > represents an atomic orbital on site (m, n). 
Indexes run from 1 to L, and the symbol <> denotes that 
the sum is restricted to the existing nearest-neighbors of 



site (m,n). Using Landau's gauge the hopping integral 
is i m ,n;m',n' = exp (2iri ^ L m l)2 for m = m', and 1 
otherwise. Therefore, the difference between our Hamil- 
tonian H and the one corresponding to an ideal L x L 
cluster on the square lattice is the absence of hopping to 
and from L sites chosen at random among the L 2 sites 
defining the lattice. A full discussion of the properties of 
this model for the case of a closed system and zero field 
can be found in Ref. [tPl. 



B. Conductance 

The conductance (measured in units of the quantum of 
conductance Go = e 2 //i) was computed by using the im- 
plementation of Kubo formula described in Ref. |^0| (ap- 
plications to mesoscopic systems can be found in Refs. 
pl[ and (2^]). For a current propagating in the re- 
direction, the static electrical conductivity is given by: 



G = -2 



Tr 



(hv x )Img(E)(hv x )ImG(E) 



(7) 



where Im Q(E) is obtained from the advanced and re- 
tarded Green functions: 
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lmg(E) = - Q (E) — G (E) 



(8) 



and the velocity (current) operator v x is related to the 
position operator x through the equation of motion hv x = 

H,x , H being the Hamiltonian. 

Numerical calculations were carried out connecting 
quantum dots to semiinfinite leads of width W in the 
range 1— L. The hopping integral inside the leads and 
between leads and dot at the contact sites is taken equal 
to that in the quantum dot (ballistic case). Assuming 
the validity of both the one-electron approximation and 
linear response, the exact form of the electric field does 
not change the value of G. An abrupt potential drop at 
one of the two junctions provides the simplest numerical 
implementation of the Kubo formula |2(| since, in this 
case, the velocity operator has finite matrix elements on 
only two adjacent layers and Green functions are just 
needed for this restricted subset of sites. Assuming this 
potential drop to occur at the left contact (Ic) side, the 
velocity operator can be explicitly written as, 

w 

ihv x =-J2(\lc,j >< - >< lc,j\) (9) 

3=1 

where (\lc, j > are the atomic orbitals at the left contact 
sites nearest neighbors to the dot. 
Green functions are given by: 



[EI-H- Ei(£) - E 2 (E)]G(E) = I 



(10) 
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where Ei^-E) are the selfenergies introduced by the two 
semiinfinite leads |^3| . The explicit form of the retarded 
selfenergy due to the mode of wavevector k y is: 

E(25) = \{e- e(k y ) - i^/i-(E-e(k y ))A , (11) 

for energies within its band \E~e(k y )\ < 2, where e(k y ) = 
2cos(k y ) is the eigenenergy of the mode k y which is quan- 
tized as k y = (n^ 7r)/(W+l), nk y being an integer from 1 
to W. The transformation from the normal modes to the 
local tight-binding basis is obtained from the amplitudes 
of the normal modes, < n\k y >= y/2/JW + l)sin(nfc y ). 
Note that in writing Eq. (|Tl]) we assumed that the mag- 
netic field was zero outside the dot p3] |, 

C. Numerical Procedures 

Input/output leads were attached at opposite corners 
of the dot as follows: input lead connected from site (1,1) 
to site (1,1 + W), and output lead from site (L, 1) to site 
(L, 1 + W). We have checked that changing the sites at 
which leads are attached does not qualitatively modify 
the results discussed here. The conductance was aver- 
aged over disorder realizations (local distribution of va- 
cancies) and within selected energy ranges. The latter 
were chosen to fit the number of channels in the leads. 
More specifically, for leads with iV c h channels energy av- 
erages were carried in the range, 

E e i E N ch ,E Nch+1 ] , (12) 

for iV c h channels in the leads, where, 

/ TTTl \ 

E « = - 2 { 1 + C0S WTT) ■ (13) 

Some calculations were also carried out at a fixed Fermi 
energy. In all cases averages were done over at least 1200 
values of the conductance. 

III. RESULTS 

A. Zero field conductance 

Fig. |l| shows relative deviations of the conductance 
with respect to the RMT result (see Eq. (Bh) for nar- 
row and rather wide leads as a function of the dot size L. 
It is noted that for small W deviations are always smaller 
than 5%, and typically below 2%. The results fluctuate 
more appreciably for the narrower lead (W=l) as ex- 
pected p2j. Relative deviations from the RMT result 
are significantly larger for W = 9 and 18. The results of 
Fig. [l] suggest that the difference with respect RMT is not 
a size effect. The larger deviation, and stronger variation 



in the explored range of L, observed for 7V c h = W = 9 is 
likely a consequence of the important contribution that 
the center of the band [E = 0) has in that case. Both 
the center of the band and its bottom [E = —4) show 
rather odd behaviors. In particular at E — no weak 
localization effect was observed (see below). 

The change in the zero field conductance as the num- 
ber of channels is varied, for fixed dot size, is illustrated 
in Fig. |[ The results for iV c h = W/2 can be accurately 
fitted by means of a straight line, as expected, although 
the slope is smaller than the RMT prediction (see caption 
of Fig. |2| and Eq. EJ) . The slope of the straight line varies 
with the ratio N^/W, or alternatively the average Fermi 
energy; for instance at E = it is actually larger than 
0.5. For fixed leads width W — 22 and a variable number 
of channels a large deviation with respect to a straight 
line is instead observed. This deviation, which is a con- 
sequence of the concomitant change in the Fermi energy 
as the number of channels is varied, increases with the 
number of channels, likely due to the increasing impor- 
tance of the contribution of the band center. Numerical 
results indicate that at the band center the conductance 
shows a much stronger dependence (increase) on the dot 
size that at any other energy within the band, probably 
due to the building up of the singularity in the density 
of sates characteristic of the square lattice at that en- 
ergy. These results suggest that if the N c ^ dependence 
has to be investigated it is more reliable to work at a 
fixed N c h/W ratio and vary the leads width. 

B. Magnetoconductance: Weak Localization 

We first discuss the energy dependence of the critical 
flux and of the weak localization term. This was done 
by investigating rather large W and varying the number 
of channels in each lead. This is equivalent to vary the 
energy range over which the energy was calculated (see 
above). Fig. || depicts numerical results for cavities of 
linear size L=78 and leads of width W= 22 (a) and 44 
(b) . The conductance was obtained by averaging over 60 
disorder realizations and 21 energies in the ranges cor- 
responding to the number of channels in the leads (see 
subsection IIC). The weak localization peak shows the 
expected behavior. The numerical results were fitted by 
means of Eq. ([!]) . At this stage it is worth noting that the 
conductance remains constant in a wide range of fluxes 
only for small W . For large W the conductance follows 
Eq. (Q) in a rather narrow range of $ and then in- 
creases steadily. This deviation from the Lorentzian-like 
law hinders the fitting of the numerical results. The fit- 
ted parameters are reported in Table I. The parameters 
derived from RMT (Eqs. (|^])) are also given in the 
Table. Eq. (||) was used with k = 1 as its explicit de- 
pendence on energy was not given in Ref. H]; this will 
suffice, however, to illustrate our point concerning the 
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strong energy dependence of that constant. We first note 
that the weak localization term is smaller than the val- 
ues predicted by RMT likely due to the large values of 
W (see below). However, the dependence of a/b on the 
number of channels is the correct one (it increases with 
N ch ) but for W = 44 and N ch = 28. The latter devi- 
ation is a consequence of the increasing importance of 
the band center. At that energy the results indicate that 
G($) decreases as a function of the flux, i.e., there is no 
weak localization effect. Due to computing limitations 
we have not been able to check whether this result is a 
size effect. The numerical results for parameter b indi- 
cate that it increases with the number of open channels, 
a behavior opposite to that given by RMT with k= con- 
stant. This suggests that it cannot be safely concluded 
that the critical flux is proportional to the square root 
of the number of open channels, as, increasing the Fermi 
energy, not only increases N c h but it also dramatically 
changes constant k in Eq. (||) . The energy dependence of 
k already appears in the supersymmetric cr-model result 
jiof . We have checked that if the energy dependent factor 
in Eq. (|^), namely, vDq, is included (the mean free path 
was calculated following the procedure of |l3|,^l| , see also 
p4|) the dependence of the RMT result on the number 
of channels is reversed, in agreement with our numerical 
results. 

The effects of the dot size on the weak localization peak 
were investigated for L in the range of L = 27-137 and 
three combinations of (W, AT c h). Averages were identi- 
cal to those mentioned in the preceding paragraph. The 
results are depicted in Figs. ^ and (|J). The weak localiza- 
tion term (a/6) shows a slight size dependence at small 
L, saturating for L approximately larger than 50 (see 
Fig. ^). This indicates that the smaller values of a/b ob- 
tained in our calculations, with respect to RMT, is prob- 
ably not a size effect. The results for (W, A r c h)=(l,l) are 
slightly smaller than those for (2,1) surely due to the con- 
tribution of the band center in the first case. The results 
for (10,5) are larger than the other two, in agreement 
with RMT. On the other hand our results for constant b 
increase with L as expected (see Eq. (0)). The numeri- 
cal results can be reasonably fitted by means of straight 
lines as shown in Fig. |^. The differences in the slopes is a 
consequence of the energy dependence discussed above. 

In order to get rid as much as possible of the strong 
energy dependence of the shape of the weak localiza- 
tion peak, we have carried out the study of the effects 
of the leads width at a fixed energy. We have chosen 
E = —2.001 (away from the band center and bottom) 
which approximately correspond to A c h = W/2. We 
fixed the dot size at L — 78 and varied the leads width 
in the range W=4-78. The results are shown in Figs. ||- 
The conductance versus the magnetic field for small 
and large W is depicted in Fig.(^|). It is readily noted 
that both the weak localization term and constant b (or 
the inverse of the square root of the critical flux) sharply 



decreases with W. Although the results are nicely fitted 
by means of Eq. |l|) the deviation of the numerical re- 
sults with respect to that equation which occurs at large 
W (see above) is already observed for W — 78 (note 
that the fitting closely follows the numerical results only 
up to $ w 1.5). The decrease of b with W, or, alter- 
natively, with jVch is illustrated in Fig. 0. The results 
can be satisfactorily fitted by means of the RMT result 
(see caption of Fig. f7|). On the other hand the weak 
localization term shows a size dependence that has not 
been previously anticipated. At small W (or number of 
channels) it increases as predicted by Eq. (||). However, 
beyond W ~ 0.2L it begins to decrease sharply reaching 
a value slightly larger than 0.05 for W = L. To explore 
the possibility that the weak localization term vanishes 
in the large L limit we have calculated the magnetocon- 
ductance for W = L, E = —2.001 and L in the range 
30-126. The numerical results were fitted by means of 
Eq. (|l|) with the parameters reported in Table II. The 
results clearly indicate that a/b does not vanish as L in- 
creases. The fact that b is almost independent of L is a 
consequence of the dependence of b on the ratio L/N^ 
(note that by taking W = L and a fixed energy the num- 
ber of channels is proportional to L). 

IV. CONCLUDING REMARKS 

Summarizing, we have presented what we believe to be 
the first detailed numerical study of the effects of Fermi 
energy, leads width and dot size on the shape of the weak 
localization peak in quantum chaotic cavities. The study 
was carried out on a model that was recently proposed 
by us |l7| which shows all the expected features of closed 
chaotic quantum billiards. Although the conclusions of 
our investigation qualitatively agree with most predic- 
tions of random matrix theory, some significant issues 
have to be highlighted. We first note that our results 
show that albeit the critical flux is proportional to the 
square root of the number of open channels, as predicted 
by RMT, the proportionality constant strongly depends 
on the Fermi energy in agreement with Efetov's analysis 
|Tof . This introduces a model (system) dependence which 
makes theoretical (experimental) comparisons with RMT 
rather delicate. Our results clearly illustrate the size de- 
pendence of the critical flux, in particular $ c cx 1/ ^/L, in 
agreement with Efetov results jl(| and the RMT results 
reported in Refs. jl5],[l6| (note that this size dependence 
was not found in a previously published RMT study, see 
Ref . [H ) . Finally, we have investigated the effects of the 
leads width concluding that the weak localization term 
sharply decreases with the ratio W/L (a result that has 
not been previously reported) , although it is likely that it 
does not vanish in the infinite L limit. This suggest that 
RMT is probably not valid for sufficiently open systems. 
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TABLE I. Fittings of numerical results such as those of 
Fig. |^ by means of Eq. (|l|) for 78 x 78 chaotic cavities with 
leads of widths W—22 and 44 attached at opposite corners 
of the dot, as discussed in the text. Chaoticity was induced 
by introducing L vacancies within the bulk. The magneto- 
conductance was obtained by averaging over 60 disorder re- 
alizations and 21 energies in the ranges corresponding to the 
number of channels N c h in the leads (see text). RMT re- 
sults as obtained from Eqs. (Q) and (|H|) with k = 1 are also 
reported. 
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5.12 
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0.151 
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0.24 
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0.98 


1.83 


0.076 


0.88 


0.22 
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4.62 


2.46 


0.097 
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8.47 


2.81 


0.125 


0.18 


0.24 




28 


12.87 


4.73 


0.091 


0.14 


0.25 



TABLE II. Fittings of numerical results such as those of 
Fig. ^| by means of Eq. (jl|) for LxL chaotic cavities with leads 
of widths W = L attached at opposite corners of the dot, as 
discussed in the text. Chaoticity was induced by introducing 
L vacancies within the bulk. The magnetoconductance was 
obtained by averaging over 1260 disorder and at a fixed energy 
£=-2.001 (which roughly corresponds to iV c h = W/2). 



L G(0) b a/b_ 

30 6.50 0.66 0.058 

42 9.34 0.89 0.051 

54 12.19 0.91 0.056 

66 15.04 0.76 0.070 

78 17.91 0.85 0.059 

102 23.65 0.98 0.053 

114 26.49 0.87 0.054 

126 29.37 0.88 0.047 
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FIG. 1. Relative deviations of the numerical results for 
the zero field conductance with respect to the RMT result 
versus dot size L. Leads of width W were attached at opposite 
corners of the dot. Chaoticity was induced by introducing L 
bulk vacancies (see text). The results correspond to averages 
over 60 disorder realizations and 21 energies in the ranges 
corresponding to the number of channels 7V c h in the leads (see 
text). The results correspond to (W,N c h) = (1,1) -diamonds-, 
(3,3) -squares-, (18,9) -circles- and (9,9) -triangles-. The lines 
are guides to the eye. 
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FIG. 3. Magnetoconductance as a function of the mag- 
netic flux (in units of their respective quanta) in 78 x 78 
chaotic cavities with leads of width W = 22 (a) and 44 (b) at- 
tached at opposite corners of the dot, as discussed in the text. 
Chaoticity was induced by introducing L — 78 bulk vacancies 
(see text). Averages were taken over 60 disorder realizations 
and 21 energies in the ranges corresponding to the number of 
channels iV c h in the leads. The numerical results correspond 
to iV c h = 4 (triangles) and 12 (squares) and were fitted by 
means of Eq. (|l|) with the parameters reported in Table I. 
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FIG. 2. Zero field conductance G(0) in units of the con- 
ductance quantum in dots of linear size L=95 versus the num- 
ber of channels in the leads iV c h- Leads of width W — 22 
(empty circles) and W = 2iV c h (filled circles) were attached 
at opposite corners of dots of linear size L = 78. Chaotic- 
ity was induced by introducing L bulk vacancies (see text). 
The results correspond to averages over 60 disorder realiza- 
tions and 21 energies in the ranges corresponding to the num- 
ber of channels iV c h in the leads (see text). The straight 
line (broken curve) fitted to the results for W = 2N C ^ is: 
G(0) =0.44iV ch -0.15. 
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FIG. 4. Weak localization term (a/b in Eq. (|lj)) as a func- 
tion of the dot size. Leads of width W were attached at 
opposite corners of the dot. Chaoticity was induced by in- 
troducing L bulk vacancies (see text). Averages were taken 
over 60 disorder realizations and 21 energies in the ranges cor- 
responding to the number of channels jV c h in the leads (see 
text). The results correspond to (W,N c h) = (10,5) -circles-, 
(2,1) -triangles- and (1,1) -squares-. The lines are guides to 
the eye. 
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FIG. 7. Constant b in Eq. (|fj) as a function of the number 
of channels in the leads . Leads of width W were attached 
at opposite corners of dots of linear size L=78. Chaoticity was 
induced by introducing L bulk vacancies (see text). The re- 
sults correspond to averages over 1260 disorder realizations, 
W=4-78 and a fixed energy _E=-2.001 (which roughly corre- 
sponds to N c h = W/2). The broken line is the RMT result 
obtained from Eq. (pi) with k — 15. 
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FIG. 5. Same as Fig. (j4J) for constant b in Eq. (hj). The 
fitted straight lines are also shown. 
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FIG. 6. Magnetoconductance as a function of the mag- 
netic flux (in units of their respective quanta) in 78 x 78 
chaotic cavities with leads of width W — 4 and 78 (squares 
and circles, respectively) attached at opposite corners of the 
dot, as discussed in the text. Chaoticity was induced by in- 
troducing L = 78 bulk vacancies (see text). The results corre- 
spond to averages over 1260 disorder realizations and a fixed 
energy _E=-2.001 which roughly correspond to N c h = W/2. 
The numerical results were fitted by means of Eq. (hi) 
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FIG. 8. Weak localization term (a/b in Eq. (|lj)) as a func- 
tion of the leads width W . The leads were attached at op- 
posite corners of dots of linear size L—78. Chaoticity was 
induced by introducing L bulk vacancies (see text). The re- 
sults correspond to averages over 1260 disorder realizations, 
M / =4-78 and a fixed energy i5=-2.001. The line is a guide to 
the eye. 
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